Light-induced current in molecular junctions: Local field and non-Markov effects 
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We consider a two-level system coupled to contacts as a model for charge pump under external 
laser pulse. The model represents a charge-transfer molecule in a junction, and is a generalization 
of previously published results [B. D. Fainberg, M. Jouravlev, and A. Nitzan. Phys. Rev. B 76, 
245329 (2007)]. Effects of local field for realistic junction geometry and non-Markov response of the 
molecule are taken into account within finite-difference time-domain (FDTD) and on-the-contour 
equation-of-motion (EOM) formulations, respectively. Our numerical simulations are compared to 
previously published results. 
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I. INTRODUCTION 

Driven transport and coherent control at the nanoscale 
are well established areas of research. Quantum 
ratchets,—!^ molecular charge,— spin^i^ and heat pumps,— 
and nano-plasmonics^ are just several examples of areas 
of recent developments. Advances in optical techniques, 
in particular near-field optical microscopy, allow single 
molecule manipulation^ and induction of bond specific 
chemistryiiS Combined with molecular junction fabrica- 
tion techniqueSfii optical spectroscopy methods are be- 
coming an important observation and diagnostic tool in 
molecular electronics J^"— 

Experimental developments led to surge of theoretical 
activity in the field of optically assisted transport^"— 
and optical response of molecular junctions.—"— 

In particular, Ref. [l^ considered molecular junctions 
composed of molecules with strong charge-transfer transi- 
tion into their excited state^^"— as a possible constituent 
for light-induced molecular charge pump, when change 
of molecular dipole occurs along the junction axis. Con- 
sideration was done within a two-level (HOMO-LUMO) 
model with ground and excited (HOMO and LUMO) 
states of the molecule strongly coupled to different con- 
tacts. In junction setup optical excitation brings elec- 
tron from occupied ground to empty excited state, and 
asymmetry in coupling to contacts assures appearance of 
current. The model was treated within non-equilibrium 
Green function approach, and perturbation theory in 
coupling to laser field was employed. 

Later Ref. [29l generalized the consideration of Ref. [l6l to 
strong laser fields. Pumping optical field was treated as a 
classical driving force, and closed set of EOMs for observ- 
ables (electronic populations and coherences of the levels 
and single time exciton correlation function) was formu- 
lated. One of the most important advances in Ref. [2^ 
was consideration of chirped laser pulses, which allowed 
formulation of charge transfer between ground and ex- 
cited states in terms of Landau-Zener problem. Chirped 
laser pulses enable to produce complete population in- 



version in molecular systems (a molecular bridge) where 
the well-known 7r-pulse excitation^ fails. 

In realistic molecular junctions optical field driving the 
molecule is a local field formed by both incident radiation 
and scattered response of the system (mostly plasmonic 
response of metallic contacts). Another feature of molec- 
ular junctions is hybridization of states of a molecule with 
those of contacts. The latter leads to non-Markov effects 
in response of the junction. 

In this paper we generalize studies reported in Ref. [2^ 
incorporating the aforementioned effects into consider- 
ation. Dynamics of local electromagnetic fields is sim- 
ulated within the FDTD technique for realistic geom- 
etry of a molecular junction similar to our previous 
publicatioui^ Non-Markov effects of junction response 
are introduced within non-equilibrium Green functions 
equation-of-motion (NEGF-EOM) approach. 

Structure of the paper is the following. After introduc- 
ing the model in section |lll we describe a junction geom- 
etry and numerical approach used in calculations of local 
electromagnetic fields in section Hill Section HVl discusses 
calculation of local field- induced electron flux through the 
junction, and section |V] introduces set of NEGF-EOMs. 
Numerical results and discussion are given in section IVll 
Section IVIII summarizes our findings. 



II. MODEL 

A model junction consists of a molecule coupled to 
two metallic contacts driven by external radiation field. 
The radiation is a time-dependent local electromagnetic 
field E{t) calculated within FDTD technique for bowtie 
geometry of the contacts (see section IIIII for details). 
Molecule is represented by a two- level system |1 > and 
1 2 > (HOMO and LUMO or ground and excited states), 
and is placed in a 'hot spot' of the local field. Contacts 
L and R are assumed to be free charge carrier reservoirs, 
each at its own equilibrium. Difference in their electro- 
chemical potentials defines bias applied to the junction 
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FIG. 1: (Color online) A sketch of local field driven molecular 
charge pump. 



eV = — fJ-R- Following Refs. Il6ll29l we consider two 
types of coupling between molecule and contacts: charge 
and energy transfer. Hamiltonian of the junction is 



Hit) =Hoit) + V 



eknk 



(1) 
(2) 
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Tn=l,2/ce{L,_R} 
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Here dj„ {dm) and c[. (cfc) are creation (annihilation) op- 
erators of electron in level m of the molecule and in state 
k in the contact(s), respectively, n,„ = a[^dm is the op- 
erator of electronic population in level to, £^12 = '^^ 



d[d2 



is operator of molecular de-excitation (£'21 = ^12)1 ^^'i 
/i is molecular transition dipole moment. Terms on the 
right-hand side of ^ represent molecular structure (two- 
level system), contacts, and coupling to the driving field. 
Right-hand side of Eq. ([3]) introduces electron and energy 
transfer between molecule and contact(s). Eqs. (H])-© 
introduce the model of Ref. [2^ with rotating wave ap- 
proximation relaxed, and with driving force treated as a 
local electromagnetic field. 

To simulate molecules with strong charge-transfer 
transition with dipole moment oriented along the junc- 
tion axis below we assume that ground state (or HOMO), 
|1 >, is coupled strongly to the left contact L, while ex- 
cited state (or LUMO) - to the right contact R. Such 
setup works as a local field driven charge pump (see 
Fig. [T|). Note that similar selective coupling can also 
be obtained for the bridge made of a quantum dots as 
discussed in Refs. ISlMssl 



III. LOCAL FIELD SIMULATIONS 

Calculations of the local electromagnetic field dynam- 
ics are carried out utilizing FDTD techniquci^ Following 
Ref. [29I we assume that the incident field, i?jj-^(,(t), has the 
form of a linear chirped pulse 



E-,j,At) = Re £0 cxp 
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— iujQt 
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where Sq is the incident peak amplitude, ujq is the inci- 
dent frequency, and parameters S and p, describing inci- 
dent chirped pulse are given by 
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with To = tpQ I \J2 log 2 (the value of the pulse duration 
ipo of the corresponding transform-limited pulse used in 
simulations is 9.34 fs) and $"(ajo) is the chirp rate in 
the frequency domain. Throughout the simulations the 
incident field is taken in the form of ([U and is normalized 
to preserve the total energy of a laser pulse at different 
chirp rates according to 
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dtEfj^^{t) = const. 



(7) 



The geometry considered is depicted in the inset of 
Fig. showing the top view of the bowtie antenna. 
To investigate the influence of chirped incident pulses on 
plasmon dynamics we choose incident field in the form 
(HI) and vary $"(i'o) = 'in'^^" {ojq). Below we shall write 
<t>". having in mind ^"(i^q). Wc further presume that 
the incident pulse is x-polarized and propagates along 
the z axis with the incident frequency at the plasmon 
resonance (sec the inset of Fig. [2^). Material dispersion 
of silver is taken in the Drude form with other numerical 
parameters as in Ref. [2^ For a given set of material and 
geometric parameters the local electric field enhancement 
exhibits well pronounced plasmon resonance as seen in 
Fig. [2^ reaching the value of 2800 near 2 eV. 

Our goal is to take plasmonic effects (local field en- 
hancement and phase accumulation) directly into ac- 
count and investigate how such crafted local fields affect 
transport properties of molecular junctions placed in the 
gap of bowtie antennas. However it is informative first 
to examine general features of chirped pulses interacting 
with plasmonic materials. It has been noted in several 
papers^^"— that local field enhancement depends sensi- 
tively on the sign of chirped excitation pulses. Moreover 
careful examination of spatiotemporal dependence of lo- 
cal fields on chirp rates^ revealed a complex dynamics 
of plasmon wavepackets that are noticeably infiuenced by 
chirped laser pulses - one may find different local points 
for a given plasmonic system, where positive chirps lead 
to higher local fields and the other way around. 
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FIG. 2: (Color online) Results of FDTD simulations for 
chirped pulses exciting the bowtie antenna schematically de- 
picted in the inset of panel (a). Panel (a) shows the enhance- 

I 

ment of the local intensity \E\ detected in the gap of the 
bowtie antenna as a function of the incident frequency. Panel 
(b) presents local field component of electric field along the 
axis of symmetry of the structure, Ex, as a function of time 
for two chirp rates at the plasmon resonance, loq = 2.057 
eV: black solid line - ^"{vq) = -3000 fs^, red dashed line - 
$"(i-'o) = 3000 fs^. Panel (c) shows the amplitude of the local 
field in the frequency domain (note that it is independent from 
the phase rate). Panel (d) represents the phase of the local 
field in the frequency domain at two chirped rates: black solid 
line - ^"{uo) = -3000 fs^ red dashed line - $"(i'o) = 3000 
fs^ 



Generally speaking, plasmonic materials can be consid- 
ered as pulse shapers^ due to high material dispersion 
near plasmon resonance, which induces a phase in the 
frequency domain resulting in shaping of the total elec- 
tromagnetic field in time domain. This is illustrated in 
Fig. [2)D-d, where one can clearly see that the positive 
chirp leads to the compression of the local field (Fig. [Jb) 
and hence stronger field enhancement. While the field 
amplitude in the frequency domain is not affected by the 
chirp sign (Fig. [2}:), obviously the phase of the field is 
significantly different for positive and negative chirp as 
shown in Fig. [2li. We note that one can not recover 
data obtained for negative chirp, for instance, by simply 
fiipping the sign of the phase for the positive chirp. Ad- 
ditional phase induced by the plasmonic system, which 
depends on the sign of the chirp rate, makes this problem 
time irreversible!^ 



IV. CURRENT THROUGH THE JUNCTION 

Time-dependent current through the junction under 
external driving is^iiS 



lKit)^-j:{Tr[T^pit)] (8) 
+ -Im y dEfK {E) Tr [r^ C- (t , i;)] j 

Here the trace is taken over molecular subspacc, is 
matrix of electronic decohcrence due to coupling to con- 
tact K 



r^„, = 27r ^ Vkm' S{E - Ek) 



(9) 
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which is energy independent in the wide-band approx- 
imation, /kIE) is Fermi-Dirac thermal distribution in 
the contacts, p(t) = —iG^{t,t) is non-eqiulibrium re- 
duced density matrix of molecular subsystem, G^'^{t,t') 
are matrices in molecular subspace of retarded and lesser 
projections of single-electron Green function 



G'mm'(r, r') = -i{Tcdm{T)dl^,{T')) 



(10) 



(Tc is contour ordering operator), and G^'{t,E) is the 
right side Fourier transform of the retarded Green func- 
tion 

/ + 00 
dt' e'^^'-'"> G''it,t') (11) 
-oo 

We are interested mostly in effectiveness of the device 
as a charge pump, i.e. we will calculate excess charge 
transferred through the system during the laser pulse 



Qif(t)- r at' {iK{t')~ii^\ 

J — oo 



(12) 



where Ixit) is defined in Eq.dS]) and if^ is current at 
bias induced steady-state condition, i.e. in the absence 
of radiation - E{t) = 0. 



V. EQUATIONS OF MOTION 

Markov approximation employed in Ref . [2^ comes from 
consideration of time-local quantities only. This ap- 
proach is sufficient when one can neglect broadening of 
molecular states induced by hybridization with states in 
the contacts. In realistic molecular junctions such hy- 
bridization is non-negligible, since molecules are usually 
chemisorbed on at least one of the contacts. Here (in 
addition to local field formation) we are going to ex- 
plore how non-Markovian effects influence characteristics 
of laser pulse induced charge pumping. 

To keep non-Markov effects a time-nonlocal quantity - 
single-particle Green function, Eq. lfTUl) - is at the focus 
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of our consideration. We employ Keldysh contour based 
EOM approach, similar to the one employed in our earlier 
publication^ (see Appendix [Xl for derivation) 

d 

i-Q^Gmm'iT,T') = 6m,m'S{T,T') + SmGmm' (t, t') (13) 

{t,ti)G (n,T') 

mi "'c 

^ Qinrh-i ,m'mi ('^i '^1 ; 7 '^l^") 



and Q is molecular subspace two-particle Green function 

{ti,T2;t3,T4) = (16) 
, Tc dm, (ri) dm^ (^2) ('^4) (^-3) , 



Here m is molecular level other than m (e.g. for m = 1 
rfi — 2), gk{T,T') is a single-particle Green function of 
free electron in the contacts 



gUr,r')^-^(T,c,{r)cliT')) 



(14) 



Smm' (t, t') = Y.K=L,R. ^mm' ('^^ is the sclf-encrgy due 
to coupling to contacts with 



^^ra'ir,r')=Y,Vrnkgk{T,T')VK 



km' 



(15) 



Note, in derivation of (|T3)) wc treated the energy transfer 
term, Eq.(l3l), at the second order of the perturbation 
theory. 

Presence of many-body interaction does not allow to 
close hierarchy of equations exactly. To make the prob- 
lem tractable we employ Markov approximation in treat- 
ing energy transfer, last term on the right in (|13p . and 
in writing EOM for two-particle GF (see below). These 
approximations are similar to those introduced previ- 
ously in Refs. I20il29l . Molecule-contact coupling in (fTS]) 
is treated exactly, thus introducing non-Markov effects 
into description. This leads to system of equations (see 
Appendix |X] for derivation) 
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I 



(17) 
(18) 

(19) 



(20) 



Here £mm = Em - em, nm(t) = pmm{t) [m ^ 1,2) are 
populations of molecular levels, p{t) = P2i{t) is molecular 

coherence, Nuit) = (^D^t)D{t)'j = e^i2,i2(t+, i, t+) is 
the molecular excitation correlation function, Tmm' = 
T1,K=L R ^mm' ^hc matrix of electronic decohcrence 
due to electron transfer between the molecule and 
contacts, with F,^„^, defined in Eq.®, T^^^'i^) = 
J2k=l R^mm''^iE) greater (lesser) projections of self- 



energy due to coupling to contacts with 



-'ram' 



[E) ^iV^,,,fK{E) 



^mm'(-^) = ^ i^-mm' [1 ^ fK{E)] 



(21) 

(22) 



and B{E) is the dissipation rate due to energy transfer 



B{E)^2n J2 



E 

K=L,Rki^k2&K 
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Note that in p7|) we omitted term coming from en- 
ergy transfer, since contribution to the total retarded 
self-energy E*" from molecule-contacts electron transfer 
~ r is much bigger than corresponding contribution 
from energy transfer ~ B{e2i)n2 (~ i3(e2i)(l — ni)) for 
m = 1 (m = 2) in a reasonable parameter range.— 
EOMs p7|) - (Pm) form a closed set of timc-depcndcnt 
equations to be solved simultaneously on energy grid 
starting from a steady-state initial condition correspond- 
ing to biased junction before the laser is switched on. 
Density matrix p{t) and retarded GF G^^t^E) obtained 
as the solution are used in ([5]) and (fT^ to calculate time- 
dependent current and excess charge pumped through 
the junction, respectively. 

In the limit of weak molecule-contact coupling F — > 0, 
neglecting local field and non-Markov effects, disregard- 
ing off-diagonal terms in spectral function, and assuming 
rotating- wave approximation Eqs. pTjl - ipU)) are reduced 
to results of Ref. [2^ (see Appendix |B] for details) . 



VI. RESULTS AND DISCUSSION 

Here we present results of numerical simulations for the 
model (HI)-® with local field formation and non-Markov 
effects taken into account as described above. Time de- 
pendent local electromagnetic field is calculated solving 
Maxwell's equations on a grid (see section llTl|) for metal- 
lic contacts of a bowtie geometry. 

Molecule is placed in a 'hot spot' situated between the 
contacts, and local field plays a role of external driv- 
ing force in electronic calculations (as described in Sec- 
tion |V|. Unless stated otherwise parameters of the elec- 
tronic simulations: temperature is 300K, molecular elec- 
tronic level positions £i = — leV and £2 = leV, elements 
of electronic decoherence matrix are T^^ = V22 — 0-leV, 
^22 = = O.OleV, T[f = F2i^ = 0, coupling to 
external field jiEo — O.OOScV (after normalization ^ for 
$" = 2000 fs^; also below coupling to external field below 
is given renormalized according to ([7]) for particular $"). 
Fermi energy is taken as origin Ep = 0, bias is applied 
symmetrically ^l,r = Ep± |e|T4d/2. AH calculations ex- 
cept Fig. [6] below are done at equilibrium, Vsd = 0. Only 
processes of energy relaxation on the molecule are taken 
into account with B{ei2) = and B{e2i) = O.leV. Time 
grid is taken from the external driving field simulations. 
Energy grid spans region from —20 to 20eV with step 
O.OOleV. Other parameters are introduced separately for 
each calculation. 

Figure [3] demonstrates pumped charge build-up dur- 
ing the laser pulse excitation. One sees that the local 
field formation leads to asymmetry in pumped charge for 
opposite chirp rates. Negatively chirped incoming field 
creates longer local pulse (see Fig. [^b), which results in 
increase in total charge pumped through the junction. 
Role of electron-hole excitations in the contacts on charge 
buildup is shown in Fig. [3^. Since processes of escape 
from LUMO into the right contact and energy relaxation 
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FIG. 3: (Color online) Time dependence of charge pumped 
through the junction Q(t), EQ. p2|) for several chirp rates. 
Shown are results for (a) $" = 2000 and -2000 fs^ without 
energy transfer B{e2i) = (solid line, red and dashed line, 
blue, respectively), and taking into account energy transfer 
term B{e2i) = 0.1 eV (dash-dotted line, red and dotted line, 
blue, respectively), and (b) $" = 10000 and -10000 fs^ for 
B{e2i) ~ 0.1 eV and nSo = 0.03 (solid line, red and dashed 
line, blue, respectively) and /ifo = 0.003 eV (dash-dotted line, 
red and dotted line, blue, respectively). 



on the molecule compete for the excited state population, 
current (and consequently pumped charge) diminish with 
increase of coupling to electron-hole excitations in the 
contacts. Fig. [Sb shows effect of intensity of incoming 
pulse on the transfered charge buildup. For higher inten- 
sity the build-up demonstrates saturation in the middle 
of the pulse. The reason for this behavior is the com- 
petition between timcscalcs related to Rabi oscillation 
induced by local field between molecular levels and elec- 
tronic escape rate from molecule into contacts (~ l/F). 
On the one hand, both negatively and positively chirped 
pulses in the middle have frequency approximately at 
resonance with HOMO-LUMO transition, uj £2 — £1, 
which is a prerequisite to effective electron transfer and 
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thus increase in pumped charge. On the other hand, at 
resonance Rabi osciUations^ at high enough intensities 
compete with electron escape rate, thus effectively block- 
ing current through the junction. Depending on parame- 
ters this may lead either to most effective charge transfer 
in the middle of the pulse (dash-dotted and dotted lines 
in panel b), or to suppression of charge transfer at this 
point (solid and dashed line in panel b). Note that the 
effect is not related to non-Markov relaxation, i.e. this 
behavior is observed also in the absence of hybridiza- 
tion between molecule and contact(s) states, and its rela- 
tion to Landau-Zener problem^ in terms of total charge 
pumped across the junction was discussed in Ref. |29|. 
Note also, that with positively chirped pulse changing 
frequency from lower to higher transfered charge buildup 
is more effective at the start of the pulse (at lower fre- 
quencies), while for negatively chirped pulse more effec- 
tive buildup takes place at the end of the pulse (compare 
solid and dashed lines in Figs. [3)3 andlSJa). Contrary to 
buildup suppression in the middle of the pulse, this ef- 
fect is due to molecule-contact hybridization. The latter 
leads to broadening of molecular levels, and effectiveness 
of HOMO-LUMO charge transfer depends (among other 
conditions) on integral of occupied states at HOMO and 
empty states at LUMO separated by frequency of inci- 
dent light / dEGfi{E)G22{E+bj). Clearly, at frequencies 
below resonance the latter is greater than at frequencies 
above it. 

Local field asymmetry relative to the sign of the chirp 
rate in the frequency domain leads to asymmetry in 
charge pumping contrary to symmetric situation pre- 
sented in Ref. |2l, as is demonstrated in Fig. [4^. One 
can see that the pumped charge is almost symmetric at 
high rates with asymmetry confined to the low rate re- 
gion. Difference between charge pumped through the 
junction at positive and negative chirp rates is shown 
in Fig. ^p. As discussed above duration of local field 
due to positively chirped incoming pulse is shorter than 
the one due to negatively chirped analog. This com- 
pression is the cause of less charge pumped through the 
system in the former case, which results in decrease in 
AQ = Q($") - (3(-$") in the region of $"(z^o) from 
to 3000 fs^. Indeed, at the very low rates frequency of the 
pulse does not change much, so the asymmetry is solely 
due to difference in pulse length. At a higher rates an ad- 
ditional factor appears: the most effective contribution 
to charge transfer takes place at a particular region of 
frequencies (at and just below resonance, as is discussed 
above). This region is passed quicker in the positively 
chirped local pulse, and in the region up to 3000 fs^ this 
results in increase of asymmetry, since negatively chirped 
pulse spends more time in its effective frequencies zone. 
Further increase of chirp rate leads to decrease and al- 
most disappearance of the asymmetry. The reason is 
decrease of ratio of the pulses difference to overall local 
pulse duration. 

Coupling to electron-hole excitations not only dimin- 
ishes pumped charge (compare solid and dashed lines in 
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FIG. 4: (Color online) Charge pumped through the junction 
during pulse vs. chirp rate: (a) total charge, i.e. integral of 
Q(t), Ea. (|12|l . over local field pulse duration, and (b) asym- 
metry in charge transfer between positively and negatively 
chirped incoming laser pulses, AQ = Q{\(^"\) — Q( — |<I>"|). 
Shown are results with {B{e2i) ~ O.leV, dashed line with tri- 
angles, blue) and without {B{e2i) ~ 0, solid line with circles, 
red) energy transfer. 



panel a), but also decreases asymmetry (panel b). The 
latter results from the fact that rate for molecular en- 
ergy relaxation (LUMO HOMO transition due to 
coupling to excitations in the contacts) is proportional 
to population in the LUMO (see discussion in Ref. . 
So for higher currents also energy relaxation will be more 
efficient, thus effectively compensating for the difference. 

Importance of non-Markov behavior for charge pump is 
demonstrated in Figure [5] Fig. [5^ shows pumped charge 
as function of level width (for two opposite choices of 
chirp rate). Increase in total charge pumped through 
the junction with increase in hybridization saturates at 
high strength of coupling between molecule and con- 
tacts. Such behavior is expected: at low hybridiza- 
tion there is only one frequency corresponding to reso- 
nance, where pumping is most effective, so only an 'in- 
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FIG. 5: (Color online) Dependence of charge pumping 
on molecule-contact states hybridization at fiEo = 0.03eV. 
Shown are (a) Charge pumped through the junction during 
pulse vs. electronic escape rate for chirp rate = 10000 fs^ 
(solid line, red) and —10000 fs^ (dashed line, blue), and (b) 
Normalized transferred charge build-up (charge normalized to 
total charge transferred during the pulse) vs. time for chirp 
rate $" = 10000 fs^ and Tfi = Tg, = O.leV (solid line, red) 
and leV (dashed line, blue). Red line in panel (b) is the same 
as solid red line in Fig. ^jp. 



stant' of chirped pulse eontributes to charge transfer. As 
molecule-contact coupling grows the condition of reso- 
nance transition becomes less and less strict. Eventually 
any frequency within the chirped pulse has roughly same 
effectiveness - this is the reason for saturation. Also, 
stronger coupling means more effective molecule-contact 
electron transfer, which competes more effectively with 
intra-molecular Rabi oscillation at resonance. This com- 
petition is demonstrated in Fig.[Sj3, where middle-of-the- 
pulse saturation (see discussion of Fig. [3]) disappears for 
stronger molecule-contact coupling. 

Finally, in Figure [H] we discuss influence of bias on 
charge pumping. Here we define optically pumped charge 
(excess charge) as a difference between charge pumped 



FIG. 6: (Color online) Dependence of charge pumping on bias, 
(a) Total excess charge pumped through the junction during 
pulse vs. bias for chirp rate — 10000 fs'^ (solid line, red) 
and —10000 fs'^ (dashed line, blue), (b) Excess charge build- 
up (normalized to total excess charge transferred during the 
pulse) vs. time at chirp rate = 10000 fs^ for V^d = -IV 
(dashed line, blue), OV (solid line, red), and 1 V (dotted line, 
black). Here Fn = FH = 0.5 eV. Red line in panel (b) is the 
same as solid red line in Fig. |3]d. 



through the junction with and without laser field. Fig.[6K 
demonstrates total excess pumped charge during laser 
pulse for opposite choices of chirp rate as function of bias. 
Application of bias has 2 effects on the pumping process: 

1. it depletes (populates) the HOMO (the LUMO) and 

2. it may block or release channels for electron trans- 
fer from LUMO to contact R. This leads to a situations 
when most effective optical pumping does not correspond 
to zero bias, rather we see shallow peak a.t V ^ IV. Ex- 
planation is related to the fact that broadened molecular 
levels are essentially a set of scattering channels with dif- 
ferent transmission probabilities: high conducting chan- 
nels are in the center of Lorcntzian, while channels in the 
sides of distribution are poor conductors. Optical pro- 
cess takes electron from an occupied ground state and 
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puts it into one of empty excited state. Effectiveness 
of the charge pump is defined by increase or decrease of 
current through the junction under optical pulse (see dis- 
cussion in Ref . [23l) . In particular, negative bias decreases 
effectiveness of the pump mostly due to blocking part of 
LUMO-i? escape routes. Positive bias opens additional 
escape routes at the tail of LUMO Lorentzian, facilitating 
increase in pump efficiency. However, additional effect of 
depleting HOMO and populating LUMO partially blocks 
optically induced HOMO-LUMO electron transfer, thus 
reducing overall effectiveness of optical pumping. Com- 
petition between the two proceses reveals itself as a shal- 
low peak at ~ IV. 

Time resolved charge buildup is presented in Fig ^p. 
Middle-of-the-pulsc saturation observed previously at 
equilibrium (solid line, same as in Fig. [Sb) is enhanced 
at negative (dashed line) and disappears at positive bias 
(dotted line). The reason is similar to competition be- 
tween Rabi frequency and escape rate discussed above. 
Indeed, with negative bias partially blocking fast escape 
route for the electron from excited state into right con- 
tact, Rabi oscillation plays an important role at quasi- 
resonant situation in the middle of the pulse. Positive 
opening additional routes makes Rabi oscillation less ef- 
fective. 



VII. CONCLUSION 

We consider a two-level (HOMO-LUMO) model for 
optically-driven molecular charge pump. Such pump 
may be realized as a junction formed by a molecule with 
strong charge-transfer transition between its ground and 
excited states. The junction is driven by both applied 
bias and laser pulse. The latter is treated as an external 
classical driving force. 

Our consideration is the generalization of previous 
study2^ which takes into account effects of local field ('hot 
spot' formation) and hybridization between the states of 
molecule and contact(s) (non-Markov effects). We formu- 
late approximate closed set of EOMs for single- and two- 
particle GFs. Electron transfer in the former is treated 
exactly. To close set of equations the latter are considered 
within Markov approximation. Our EOMs are reduced 
to set of equation derived in Ref. under several simpli- 
fying assumptions: weak molecule-contact coupling (ne- 
glect of hybridization), neglect of non-diagonal terms in 
molecular spectral function, and within rotating wave ap- 
proximation. 

Incoming laser pulse is assumed to be linearly chirped. 
Local field is calculated within FDTD technique on a grid 
with bowtie antenna geometry used to represent junction 
metallic contacts. We find that contrary to symmetric 
behavior of the pump relative to sign of the chirp rate, 
duration of the corresponding local field pulse depends 
on the sign of incoming chirp, which results in asymmet- 
ric operation of the pump. The asymmetry depends on 
the incoming pulse chirp rate in a non-monotonic man- 



ner. Junction response to optical driving is symmetric to 
both low and high chirp rates, going through a maximum 
between the two extremes. We find that this behavior is 
caused by the correspondence between pulse duration of 
the local field and detuning of its frequency at the end of 
the pulse from energy difference between molecular states 
(£2 - ei)- 

We note that at quasi-resonance charge pump be- 
comes ineffective, due to the competition between intra- 
molecular Rabi oscillation induced by the pulse with elec- 
tron transfer from molecule to contact. Increase of the 
molecule-contact coupling strength increases electron es- 
cape rate, thus reducing ineffectiveness of the pump due 
to Rabi oscillations. 

Also we study the effect of bias on optically-facilitated 
charge transfer through the junction. We find that in the 
non-Markov situation (i.e. when hybridization between 
molecule and contacts is non-negligible) most effective 
charge pump regime is at finite positive bias, rather than 
at equilibrium as one may expect from Markov consider- 
ation of Ref. [2^ The effect comes from optically-assisted 
charge redistribution between low and high conducting 
scattering channels in broadened molecular states, as 
was discussed in our previous publication.— Within the 
model negative bias reduces (positive increases) excess 
charge pumping due to blocking (facilitating) outgoing 
scattering channels in the excited molecular state and 
thus increasing (decreasing) the role of intra-molecular 
Rabi oscillations. 

Finally, direct electron-hole excitation in contacts, 
heating, and inelastic effects are examples of effects be- 
yond current consideration which may also have a sig- 
nificant impact on the properties of a molecular charge 
pump. 
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Appendix A: Derivation of Ea.(|13p 

EOM for ([T5)) in contour variable r starts from writing 
Heisenberg equation for dmij) 

, d , 

i-^^Gmm' {TtT ) = S,n „i'S{t, T ) + S,nGmm' {t, T ) 
OT 

- fiE{t)G^m' (r, r') + V^^Gu^, (r, r') (Al) 
+ X! ^fcTfe2^™fc2,m'fci(T-,r;r',T4-) 

The last two terms on the right come from electron and 
energy transfer terms in Eq.([3]). Treating the two within 
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non-crossing approximation allows to find the first ex- 
actly within a standard procedure^ 

(A2) 

Two particle Green function in the second term is treated 
(still keeping non-crossing approximation in mind) within 
first order perturbation theory in energy transfer 



(A3) 



dngk., (r, ri 

^9ki (j'l 1 '^)^mmi ,m'mi 

{t,ti;t',ti+) 



Using (|A4p and similar expressions for other parts of the 
Keldysh contour deformed in accordance with Langreth 
rules^ in the energy transfer term of lesser projection of 
diagonal element of (|13p . and utilizing ([^5]) leads to (|18p . 



Finally, Eq. pO|) is treated with Markov approximation 
(see Eq. (jA4p above) applied to both electron and energy 
transfer terms. Then the derivation goes along the lines 
presented in Ref l29l. 



Substituting and ((X3|) into ((XT|) yields p^ . 

Eg. lfTT)) is retarded projection of Ea. (fT3l) with omit- 
ted energy transfer term. The approximation is based 
on an estimate that in usual situation electron escape 
rate should be much bigger than corresponding energy 
transfer, T > 

Eqs. (fTSll and (fT9|) is lesser projection of (fTS]) taken 
at equal times, — iG<(t,t). Note that (fT9|) is exact, 
while in (|18p we employ Markov approximation in deriva- 
tion of the energy transfer term similar to previous 
publicationS)2£i2^ for example 

X (d\^{t)drn{t)d'l^^ (ti)dmi (ii)^ 



-OO 

X (dl(t)rf«(Ot^li(i)rf^i(i)) 

[1 - "/cal^-fciTTlS (£/c2 - £/ci + £mi ~ £mi) 
X (d7„(t)d;,,(t)dl^(t)d^,(t)^ 



Appendix B: Markov limit of Eqs. (fT7)) -(f20 ]) 



EOMs derived in Ref. |29| are Markov limit of Eqs. (fT8|) - 
(|20p within static quasiparticle approximation assumed 
for molecular states. The latter implies disregarding 
Ea. (fT7|) . and assuming instead 



i [Gl,^.{t,E) - Gl,^,{t,E)] = 2Trdm,n'S{E - £,„) (Bl) 
Then, disregarding level mixing due to coupling to con- 
tacts rf^ = = 0, Eqs. jrp and (HID reduce to 
Eqs. (33) and (34)i^ of Ref. IM After omitting non- 
diagonal elements of self-energy in Ea. ([20|) one gets 
Eq.(35) of Ref. m 
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